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ABSTRACT 

In order to better understand the possibility of coronal heating by MHD waves, we 
analyze Fe xu 195.12 A data observed with EUV Imaging Spectrometer (EIS) onboard 
Hinode. We performed a Fourier analysis of EUV intensity and Doppler velocity 
time series data in the active region corona. Notable intensity and Doppler velocity 
oscillations were found for two moss regions out of the five studied, while only small 
oscillations were found for five apexes of loops. The amplitudes of the oscillations 
were 0.4 - 5.7 % for intensity and 0.2 - 1.2 km s _1 for Doppler velocity. In addition, 
oscillations of only Doppler velocity were seen relatively less often in the data. We 
compared the amplitudes of intensity and those of Doppler velocity in order to identify 
MHD wave modes, and calculated the phase delays between Fourier components of 
intensity and those of Doppler velocity. The results are interpreted in terms of MHD 
waves as follows: (1) few kink modes or torsional Alfven mode waves were seen 
in both moss regions and the apexes of loops; (2) upwardly propagating and standing 
slow mode waves were found in moss regions; and (3) consistent with previous studies, 
estimated values of energy flux of the waves were several orders of magnitude lower 
than that required for heating active regions. 

Subject headings: Sun: corona — Sun: oscillations 
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1. Introduction 



There is much disagreement over which mechanisms heat the solar atmosphere above the 
chromosphere up to the coronal temperature of several million Kelvin. At present, two major 
possibilities are often suggested: Alfven wave heating and nanoflare heating. However, waves 
with sufficient energy flux to heat the corona have not yet been detected in the corona, except for 
waves triggered by flares. Furthermore, the typical spatial scale of nanoflares is often considered 
to be smaller than the resolution of the current telescopes. It is also uncertain which mechanism 
is dominant, or whether only one mec hanism is dominant . Although a large number of wave and 



oscillation events have been observed (lAschwan den 



20041) . to further understand coronal heating, 



it is necessary to analyze the statistical properties of waves in the corona, such as the relation 
between oscillation phenomena and their locations, and between MHD modes and their locations. 



Aschwanden et al 



(|2002|) conducted a statistical analysis of the physical parameters of transverse 
loop oscillations observed with the Transition Region And Coronal Explorer (TRACE). They 
investigated 26 coronal loops with periods of 2 - 1 1 min (decay time: 3-21 min) and concluded 
that the oscillations could be described not by the simple kink eigen mode oscillations but by 
flare-induced impulsively generated MHD waves. In another study, Doppler shift oscillations 
in hot coronal loops were investigated using the Solar Ultraviolet Measurement of Emitted 



Radiation (SUMER), in which 54 oscillations in 27 flare-like events were analyzed (|Wang et al 



2003|) . The periods of oscillations were 7-31 min (decay time: 6-37 min), which were 
somewhat longer than the TRACE loop oscillations. The authors interpreted these waves in 
terms of slow magnetoacoustic standing waves in hot loops because the phase speeds roughly 
agreed with the sound speed and the phase delays of intensity between Doppler shift were 1/4 
period. Observations of coronal oscillatory phenomena were also made with the EUV Imaging 
Spectrometer (EIS) onboard Hinode, w ith the data suggesting the existence of slow mode 



(Mar iska et al. 



2008 



Wang et al 



2008), and fast kink mode (|Van Doors selaere et al 



200 9allb|), slow mode an d fast sausage mode (|Erdelyi & Taroyan 



2008). In addition, recent studies statistically 
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showed the e xistence of propagating slow magnetoac oustic wa ves in coronal ho. 



and EIS data (IBanerjee et al 



2009; 



Gupta et al. 



2009). Finally, 



O'Shea&Dovle 



es usin g SUMER 



(12009f) conducted 



a statistical analysis of intensity oscillations in an active region using consecutive 40" slot images 
of an Fe xn 195.12 A emission line observed with Hinode/ElS. They found that oscillations with 
frequencies in the range of 2 - 8 mHz were concentrated in bright plages, while oscillations of 
higher frequencies were concentrated in the edges of those areas. The oscillations were interpreted 
in terms of magnetoacoustic waves at frequencies less than 8 mHz and fast magnetoacoustic 
waves at frequencies above that. However, their data did not have information about line profiles, 
so the identification of MHD waves was not conclusive. 

In order to reveal the nature of MHD waves in the corona, we need to take into account 
the properties of coronal structures such as loops or moss. Recent observations of active region 
loops with EIS have revealed that coronal loops are not isothermal, and that their filling factors 



are approximately 



(Warren et al 



0%, which supports the existence of multi-strands of sub-resolution size 



2008J). Moss has a bright reticulated pattern, whic h emits EUV by thermal 



conduction from overlying hot and relatively high-pressure loops (jBerger et al 



1999). The time 



scale of moss dynamics is 10 - 30 s, which is due to internal intens ity variation, obscura tion 



by chromospheric jets , or m echanical interaction with those jets (|De Pontieu et al. 



1999) 



Katsukawa & Tsunetal ((2005) found that filling factors are lower at the footpoints of hot loops 
than at the footpoints of cool loops. They concluded that braiding of the coronal magnetic field 
is more efficient at hot loops than at cool loops, which means there is a larger magnetic energy 
release at loops with a lower filling factor. Recently, an Fe xu 195 A emission line observed 
with HinodefElS was analyzed to investigate the time var iability of intensity, Doppler velocity, 



2009). There was less than 15% variability 



and nonthermal velocity in moss (|Brooks & Warrenl l: 
in intensity, 3kms 1 variability in Doppler velocity, and the average nonthermal velocity was 
21.7kms _1 . 
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In this study, we analyzed time series data of intensity and Doppler v elocity for loops and 



moss in an active region using sit-and-stare mode data observed with EIS (|Culhane et al 



2007) 



Our purpose is to identify the modes of detected oscillations in this active region in order to 
understand which modes of MHD waves contribute to the heating of active regions. We compared 
Fourier amplitudes of intensity and Doppler velocity to determine the modes of the detected 
oscillatory signals, and then investigated the correspondence between these coronal structures and 
the modes. 



2. Data reduction and analysis 

We used sit-and-stare mode data observed with EIS. The duration of the data was 2006 Dec 3 
12:43:10-13:37:20 UT, and the field of view was 1" in the ^-direction and 256" in the y-direction 
in heliocentric coordinates. There were 288 data points in the time direction. The exposure time 
was 10 s, and the average time cadence was 11.3 s, with a range of 11.3 - 11.5 s. The standard 
Solar Software interactive data language (SSWIDL) procedure eis_prep was used to calibrate the 
data. We used an Fe xn 195.12 A emission line, whose line profile at each pixel was fitted using 
a single Gaussian function. Then intensity (integrated line profile) and Doppler shift (difference 
between line centroid and 195.12 A) were calculated using the fitted parameters. After removing 
the influence of the slit-tilt and the orbital variation using the eis_wave_corr procedure, we were 
able to investigate the average line centroid at each exposure. However, a 90 - 100 min period 
component still remained, which was subtracted for accuracy. Usually, the zero-point of Doppler 
velocity is removed through examination of a quiet sun region, where Doppler velocity is much 
smaller and less dynamic. However, with a short exposure time (10 s) in our data, the incident 
photons at the quiet sun region were too small for analysis, so the zero-point of Doppler velocity 
was calculated by using all pixels along the slit. Since no explosive events were observed in our 
data, this zero-point adjustment does not produce artificial oscillations. Note that there were no 
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prominent wing c omponents in the line profiles, and the weak line blend with Fe xn 195.18 A 



(|Young et al 



2009) can be ignored bec ause the main focus of our research is the time variability 



of Doppler velocity (|Wang et al. 



2009b). 



Panel (a) in Figured] shows a 171 A image of NOAA active region 10926 taken by TRACE. 
Panel (b) in Figure Q] shows an EIS Fe xv 284.16 A raster scan image. The white vertical line 
in each panel indicates the cotemporal position of the EIS slit in both panels. In the TRACE 
image, we can see that the slit cuts across some apexes of coronal loops and some moss regions. 
Fe i 6302 A Stokes-V data from Hinode/SOT are also shown in panels(a) and (b). White and 
yellow/blue contours in the two images indicate positive and negative magnetic fields, respectively. 
SOT's field of view is represented by a white dotted rectangular area. Panel (c) in Figure \T\ shows 
the time averaged intensity at each pixel. We analyzed the structures which displayed a high 
intensity value (i.e., a high signal-to-noise ratio), and considering the correspondence between 
the structures seen in images (a) and (b) and the intensity profile of the EIS slit in panel (c), we 
identified the structures of the studied locations. A1-A5 and M1-M5 in panel (c) represent apexes 
of loops and moss regions, respectively. A small number of pixels in the top and bottom of the slit 
were excluded from the analysis due to the effect of a jitter-removing process, as outlined below. 

To remove the influence of a jitter of the spacecraft, we used data from the Ultra Fine Sun 
Sensor (UFSS). After subtracting the jitter in the N-S direction from the original time series data, 
a t rend with a time scale of ~10 min in the y-direction still remained, similar to that observed 



by 



Utzetal. 



(|2010|) . Then we calculated the cross -correlation coefficient between one shot 
and its following shot in the original data and obtained the shift in the y-direction where the 
coefficient was maximized. The absolute values of the shift from the cross-correlation method 
were somewhat larger than those of the UFSS output in the N-S direction. We shifted the original 
data using both shift values, and concluded that the cross-correlation method was more suitable. 
Using this shift process, we could remove apparent oscillations in the y-direction caused by 
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the jitter. The UFSS output in the E-W direction (x-direction) showed a maximum drift of 2". 
Bec ause this is comparable to the spatial resolution of the EIS 1" slit, which approximately equals 



2" (IHara 



20081) . we neglected the influence of the jitter in the x-direction. The processed data are 



shown in panels (d), (e), and (f) in Figured] 

Time series data and Fourier spectra for intensity, Doppler velocity, and line width at 
y = -147" (M5) in the slit are shown in Figure [21 The time series data for each pixel was 
interpolated on to a constant 11.3 s cadence grid, before being subjected to rigorous Fourier 
transforms. Panels (a), (b), and (c) show the time series of intensity, Doppler velocity, and line 
width, respectively. Before calculating the Fourier transform, we subtracted a 48-point running 
average from the original time series data in order to reduce the influence of trends longer than, 
or as long as, the duration of the data, which resulted in a decrease in power for periods longer 
than 48 x 1 1.3 - 540 s. This introduced an artificial low-frequency cut-off of 1.8 mHz, resulting in 
detections below this value being less reliable. The time series data after the subtraction process 
are shown in panels (d), (e), and (f). We defined oscillations as signals which exceed the 99% 
significance level calculated by assuming the Poisson noise of incident photons. This is indicated 
by a red line in panels (g), (h), and (i). 



3. Results 

In Figure [H pairs of Fourier components with the same frequency (I(v)/Io,v(v)/cs) are 
plotted, with the horizontal axis in each plot indicating intensity amplitude (/) normalized 
by the time averaged intensity (7 ) at each pixel, and the vertical axis in each plot indicating 
Doppler velocity amplitude (v) normalized by the sound speed (cs ~ 167 km s -1 ) at the formation 
temperature of Fe xn (~ 1.3 MK). The symbols in each plot represent data points where: only 
intensity amplitude is significant (red dots), only Doppler velocity amplitude is significant (blue 
crosses), both intensity and Doppler velocity amplitude are significant (green diamonds), and 
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neither intensity nor Doppler velocity amplitude are significant (black dots). Red triangles 
represent data points where both intensity and line width amplitudes exceed the significance level, 
while red dots indicate data points where line width amplitude is not significant. The three orange 
dotted lines show the relationship between intensity amplitude and Doppler velocity amplitude 
for magneto-acoustic waves when the angle between the magnetic field and the line of sight is 0, 
60, and 75 degrees. From Figured we found that there were only a few significant oscillations at 
the apexes of loops, while there were a substantial amount of significant intensity and Doppler 
velocity oscillations at the two moss regions. 

According to MHD wave theory, there are four major modes of magnetic tu bes in the 



corona: torsional A lfven mode, fast kink mode, fast sausage mode, and slow mode ([Roberts et al 



1984 



Spruit 



1982|). Torsional Alfven mode is incompressible in the first order, so it produces 
no density (i.e., intensity) oscillations. Assuming our observation resolves both the spatial scale 
of the torsional Alfven mode and the temporal period of the mode, then out-of-phase Doppler 
velocity oscillations can be observed at opposite edges of the waveguide. However, if either the 
spatial scale or the temporal period of the mode are unresolved, emission line broadening will 
be observed instead of Doppler velocity oscillations. Transversal oscillations of loops occur in 
fast kink mode, which means only Doppler velocity oscillations are likely to be observed as a 
marker. In addition, intensity oscillations are not observed because of the incompressibility of fast 
kink mode in the first order. In fast sausage mode, radial oscillations of loops produce intensity 
oscillations. Since the two surfaces of a loop along the line of sight oscillate in opposite directions 
and cancel each other out, emission line broadening is observed instead of Doppler velocity 
oscillations. Slow mode produces both intensity and Doppler velocity oscillations. Since the solar 
corona is an optically thin and collisionally excited atmosphere, intensity is proportional to the 
square of density (/ oc p 2 ). By performing a linear approximation, we obtained the relationship 
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between density perturbation p\ and intensity perturbation I\, 

'i = A (i) 

where /o is intensity and p is density in the unperturbed state. Then combining Eq.® with the 
Fourier transform of the linearized continuity equation, we obtained 

2/o c ph 

where Vi is velocity perturbation, and phase speed c ph is defined as c ph = co/k, where to is 
frequency and k is wavenumber. For the slow mode (c P h = c$), this relationship is indicated by the 
orange dotted line (8 = 0) in each plot of Figure |3] 

Considering the properties of the four MHD modes discussed above, we can interpret Figure 
[3] in a straightforward manner (see Tabled]). Red triangles (significant oscillations in intensity 
and line width) can be interpreted as fast sausage mode, and blue crosses (significant oscillations 
only in Doppler velocity) can be interpreted as torsional Alfven mode or fast kink mode. Green 
diamonds (significant oscillations in intensity and Doppler velocity) correspond to slow mode, and 
red dots (significant oscillations only in intensity) may be interpreted as fast sausage mode when 
line width oscillations are obscured by turbulent plasma motion, or slow mode when Doppler 
velocity signals are reduced due to the inclination of the magnetic field to the line of sight. 
Finally, black dots represent non- significant oscillations, oscillations of a shorter period than the 
time cadence, and oscillations of sub-resolution spatial scale. The above interpretations therefore 
suggest that fast sausage mode and slow mode were observed at the moss regions, while fewer 
oscillations were observed at the apexes of loops. Distinguishing between fast kink mode and 
torsional Alfven mode is possible by analyzing the phase difference at opposite edges of coronal 
structures, but this requires spatial resolution high enough to resolve the elemental magnetic 
strands. 

Phase delays between intensity and Doppler velocity were calculated for data points where 
both intensity and Doppler velocity amplitudes exceeded the significance level. The results are 
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shown in Figure HI Panel (a) shows frequency in the horizontal axis and phase delay in the vertical 
axis. Panel (b) shows the position from the bottom of the slit in the horizontal axis and phase 
delay in the vertical axis. Both panels suggest the existence of upwardly propagating slow mode 
waves (phase delay: +n) and standing waves (phase delay: ±n/2) with frequencies of 1 - 4 mHz. 
Adiabatic theory suggests that phase delays of or ±n between intensity and Doppler velocity 
oscillations are ind icative of propagating waves, while phase delays of ±n/2 suggest the presence 



of standing waves (|Deubner et al. 



1990J). Phase delays between intensity and Doppler velocity 



(i.e., V-I phase) are established for all slow mode waves. 



4. Discussion 



Assuming the coronal loops are resolved by our observations, the absence of waves at the 
apexes of cool loops means that waves excited by the photospheric granulation or nanoflares 
must have dissipated before they reached the apexes of cool loops. For the cases in which 
the spatial scales of the waves were sub-resolution size, incoherency of oscillations may have 
obscured our observat ions. This is plausi ble when considering recent observations which support 



a multi- strand model (Warren et al 



2008). In such cases, non-thermal broadening is enhanced 
in the transverse direction for torsional A lfven mode, fast kink m ode, fast sausage mode, and in 



the longitudinal direction for slow mode (Ha ra & Ichimoto 



1999). In contrast to the observed 



five apexes, our findings indicate that two moss regions out of the five studied contained large 
amounts of fast sausage mode and slow mode. Our study revealed intensity oscillations of 
0.4 - 5.7% and Doppler velocity oscillations of 0.2 - 1.2 km s _1 , with average amplitudes o f 1.4% 



and 0.5 km s , respectively. This is consistent with the results of Br ooks & Warrenl (120091) . Also 



detected were upwardly propagating slow mode waves (phase delay: n or -n) with frequencies of 
1 - 4 mHz. One possible mechanism that excites slow mode is the leakage of 5-min oscillations 



(~ 3 mHz) of the photosphere, as was reported in a study by 



De Moortel et al 



feOOCh . in which 
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they detected 2 - 4% oscillations in EUV intensity. Anoth er possible mechanism i s the mechanical 



2003). If this is the 



interaction between spicules and the surrounding corona (|De Pontieu et al. 
case, frequencies within the 3-5 mHz range might reveal certain properties of spicules, such as 
lifetime. 

We crudely estimated the energy flux at each pixel based on the interpretations in terms of 
MHD waves outlined in Section 3. After combining a linearized continuity equation and Eq. CD, 
we obtained 

to 1 1 
k 2I 

The line-of-sight effect reduces the observed Doppler velocity v obs as v obs = v\ cos 6, where 6 is the 
angle between the direction of the plasma motion and the line of sight. Using Eq. ©, the energy 
flux of compressible mode is estimated as follows: 

Po[/~i/(2/ )]V or pol^z) c s (slow mode), (4) 

\COS a I 

p [Ii/(2I )] 2 c A 3 (fast sausage mode). (5) 

The observed Doppler velocity for fast sausage mode is cancelled along the line of sight at the 
boundaries of each loop, so we used / instead of v from Eq. ©. The energy flux transported by 
torsional Alfven mode or fast kink mode is represented as 

Po(yi) 2 c A (torsional Alfven or fast kink mode). (6) 

The calculated energy flux as a function of position along the slit is shown in Figure [51 Red pluses 
and blue crosses represent fast sausage mode and fast kink mode, respectively. Green and orange 
diamonds respectively represent upwardly and downwardly propagating slow mode estimated 
using intensity amplitude. As outlined in section 3, we determined propagation direction based 
on phase delays between intensity and Doppler velocity as follows: upwardly propagating was 
defined as a phase delay in the range of (5/6)tt to n and —n to -(5/6)n, while downwardly 
propagating was defined as a phase delay in the range of — (1/6)tt to (l/6)n. Green and orange 
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triangles respectively represent upwardly and downwardly propagating slow mode estimated 
using Doppler velocity amplitude. We assumed p = 8.4 x 10~ 16 gcm~ 3 , c s = 167 km s _1 , and 
c A = 1.0 x 10 3 km s _1 . If our interpretations in terms of MHD waves are correct, our findings 
indicate that fast sausage mode waves have a significant proportion of the energy flux needed for 



Withbroe & Noyes 



( 1977 ) 



coronal heating in active regions (10 7 erg cnT 2 s -1 ), as estimated by 
For energy calculations, 6 is assumed to be zero. As a result, derived energies act as lower limits 
and rise significantly as the inclination angle of the magnetic field increases. Note that when 
Doppler velocity signals are obscured by the inclination of the magnetic field, red dots in Figure 
[3] should be interpreted as a signature of slow mode. This means that the estimated energy flux, 
shown by red pluses in Figure [5l is reduced by 2 - 3 orders of magnitude as a result of wave 
energy flux being proportional to the sound speed rather than the Alfven speed. Given that all 
red dots are below the orange line for 9 = in Figure [3l a large portion of the red dots may be 
interpreted in terms of slow mode. 



2009) 



Recent studies have detected torsional Alfven mode in the photosphere (IJess et al. 
by investigating line width oscillations of ~ 2.6 km s -1 , fast kink mode in the chromosphere 



(He et al. 



2009) with velocity amplitudes of 4.7 - 20.8 km s , and fast kink mode in the corona 



( Tomczvk et al. 



2007) with velocity amplitudes ~ 0.3 km s . Moreover, considering the extremely 



slow rate of dissipation of torsional Alfven mode and fast kink mode waves, we expected to 
detect them in our observations. However, there were few signatures of either. There are three 
possible explanations for this. Firstly, torsional Alfven mode and fast kink mode waves may have 
been reflected below the corona due to sudden changes in plasma density at the transition region. 
Secondly, they may have dissipated below the corona. However, this explanation seems to be 
less likely given their slow rate of dissipation due to their incompressibility. Thirdly, either the 
spatial scale or periods of these mode waves could have been smaller (shorter) than the resolution 
of our observations, in which case, broadening of emission line profiles occurs. To investigate 
torsional Alfven mode and fast kink mode more thoroughly, it is necessary to analyze in more 
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detail the time variability and spatial distribution of emission line width. However, this is beyond 
the scope of this paper. Such analysis will become possible with much higher spatial resolution. 
Nonetheless, we feel that our method to identify MHD wave modes provides a useful diagnostic 
tool for the statistical study of waves in the solar corona. 

Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with 
NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Scientific operation 
of the Hinode mission is conducted by the Hinode science team organized at ISAS/JAXA. 
This team mainly consists of scientists from institutes in the partner countries. Support for the 
post-launch operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, ESA, 
and NSC (Norway). We greatly appreciate the proofreading/editing assistance from the GCOE 
program. 
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Table 1: The correspondence between MHD mode waves and oscillations of observable quantities. 
Mode Density Doppler velocity Line width 

Torsional Alfven mode ^ ^ 

Fast kink mode ^ 

Fast sausage mode ^ ^ 

Slow mode £ ^0 
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(a) TRACE 171 + SOT FGIV (contour) (b) EIS 284 + SOT FGIV (contour) (c) Time averaged intensity 




12:50 13:00 13:10 13:20 13:30 12:50 13:00 13:10 13:20 13:30 12:50 13:00 13:10 13:20 13:30 

Start Time (03-Dec-06 12:43:04) Stort Time (03-Dec-06 12:43:04) Start Time (03-Dec-06 12:43:04) 



Fig. 1.— (a) TRACE 171 A image obtained at 13:36:46UT on 2006 Dec 3. (b) EIS Fe xv 284.16 A 
raster scan image obtained at 15:32: 19-17:46:3 1 UT on 2006 Dec 3 . In panels (a) and (b) , the white 
vertical line indicates the cotemporal position of the EIS slit, and Fe i 6302 A Stokes-V data from 
Hinode/SOT are also shown. White and yellow/blue contours in the two images indicate positive 
and negative magnetic fields, respectively, and SOT's field of view is shown with white dashed 
lines, (c) Time averaged intensity at each pixel on the EIS slit. The lower panels show time series 
data of the EIS sit-and-stare mode observation starting at 12:43: 10UT on 2006 Dec 3. The time 
cadence for each exposure is 1 1.4 s. (d) Intensity, (e) Doppler velocity, (f) Line width. 
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(a) Time series 



(b) Time series 



(c) Time series 
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(d) Running average subtracted 
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(g) Spectrum 
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(e) Running average subtracted 
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(h) Spectrum 
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(f) Running average subtracted 
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Fig. 2. — Examples of time series and Fourier spectra (at y = -147" in the EIS slit; M5). (a), (b), 
and (c): Time series of intensity, Doppler velocity, and line width, respectively, (d), (e), and (f): 
Residuals after subtracting a 48-point running average from the original time series data, (g), (h), 
and (i): Fourier spectra. The vertical dotted line in each plot indicates an artificial low-frequency 
cut-off of 1.8 mHz (540 s). The horizontal red line in each plot indicates the 99% significance level 
calculated from Poisson noise of incident photons. 
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Fig. 3. — Panels in the left column (A1-A5) show scatter plots at the apexes of TRACE loops, 
and panels in the right column (M1-M5) show scatter plots at the moss regions. The horizontal 
axis in each plot indicates intensity amplitude normalized by the time averaged intensity. The 
vertical axis in each plot indicates Doppler velocity amplitude normalized by the sound speed 
(cs ~ 167 km s" 1 ) at the formation temperature of Fe xn (~ 1.3MK). Red dots, red triangles, blue 
crosses, green diamonds, and black dots respectively represent data points where only intensity 
amplitude was significant, both intensity and line width amplitude were significant, only Doppler 
velocity amplitude was significant, both intensity and Doppler velocity Fourier amplitude were 
significant, and neither intensity nor Doppler velocity amplitude were significant. Three orange 
dotted lines show the relationship between intensity amplitude and Doppler velocity amplitude for 
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Fig. 4. — Phase delay between intensity and Doppler velocity at moss regions M3, M4, and M5. (a) 
Phase delay as a function of frequency. The vertical dotted line indicates an artificial low-frequency 
cut-off of 1.8 mHz (540 s). (b) Phase delay as a function of solar v. 
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Fig. 5. — Estimated energy flux in terms of MHD modes. Red pluses, blue crosses, green di- 
amonds (green triangles), and orange diamonds (orange triangles) respectively represent energy 
flux interpreted in terms of fast sausage mode, fast kink mode, upwardly propagating slow mode 
calculated by using amplitude of intensity (Doppler velocity), and downwardly propagating slow 
mode calculated by using amplitude of intensity (Doppler velocity). 



